X-Ray Flux Decrease of the Accretion-Powered Millisecond 
Pulsar SAX J1808.4-3658 in Quiescence detected by ASCA 
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ABSTRACT 

The accretion-powered millisecond X-ray pulsar, SAX J1808. 4-3658, was observed in quies- 
cence with ASCA in September 1999. We detected a dim X-ray source in the SIS data at the 
position consistent to SAX J1808. 4-3658. The source count rate was (1.1 ± 0.4) x 10"'^ cts s~^ 
(0.5-5 keV) for a single SIS, which corresponds to (3 ± 1) x 10~^* ergs s~^ cm~^ if a power 
law energy spectrum of photon index 2 with low-energy absorption corresponding to a hydrogen 
column density of 1.3 x 10^^ cm~^ is assumed. The statistical quality of the data was insufficient 
to constrain the energy spectrum or to detect the 401 Hz coherent pulsation. We compare the 
data with the BeppoSAX observation also made during the quiescent state, and find that the 
X-ray flux measured by ASCA is at least a factor of 4 smaller than that measured by BeppoSAX. 
We discuss the possible X-ray emission mechanisms that could explain the flux change, including 
the radio pulsar and the radio pulsar shock emission. 



Subject headings: Pulsars: general 
general — X-rays: stars 

Introduction 



■ Stars: individual (SAX J1808.4-3658) — Stars: neutron — X-rays: 



SAX J1808. 4-3658 is a unique transient X-ray 
source which displays both type I X-ray bursts 
and millisecond X-ray pulsations. The source 
is thought to be the formerly missing link be- 
tween the neutron-star low-mass X-ray binaries 
and the millisecond radio pulsars. SAX J1808.4- 
3658 was discovered in September 1996 by the 
Wide Field Camera aboard BeppoSAX (in 't Zand 
et al. 1998). It reached a peak intensity of 0.1 Crab 
and lasted between 6 and 40 days. Detection of 
two bright type I X-ray bursts demonstrated that 
the compact object is a neutron star. Recently, a 
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third burst was found approximately 30 days after 
the first two bursts through re-analysis of a more 
extended set of data. An improved distance es- 
timate resulted in 2.5±0.1 kpc (in 't Zand et al. 
2000). 

The source was detected again in April 1998 
with the PCA aboard RXTE (MarshaU 1998). 
The observation led to the discovery of a coherent 
X-ray pulsation at 401 Hz (Wijnands & van der 
Klis 1998). The pulse profile was sinusoidal and 
had a low fractional rms amplitude of ~4 % (2-60 
keV). Because mass accretion on to the neutron 
star was considered to power the X-ray emission, a 
major part of the accreted mass may be channeled 
to the magnetic poles without being expelled by 
the so-called propeller effect. This sets an upper 



limit to the surface magnetic field of the neutron 
star as (2 — 6) x 10^ G, comparable to that of the 
millisecond radio pulsars (Wijnands & van der Klis 
1998). Subsequent analysis showed the presence of 
a binary period of 2.0 hr (Chakrabarty & Morgan 
1998). This confirmed the low-mass nature of the 
companion. 

Several transient X-ray binaries have been de- 
tected in quiescence, and their X-ray properties 
were found to be quite different, depending on 
the nature of the compact object in the system. 
Soft X-ray transients containing a neutron star 
(NS SXTs) generally have a bottom luminosity 
around 10^^ — 10^"^ ergs s^^ in quiescence, whereas 
those containing a black hole (BH SXTs) can be 
as dim as or even dimmer than 10'^^ ergs s^^ (Asai 
et al. 1998). The former has a characteristic en- 
ergy spectrum consisting of a soft component and 
a hard tail. Quiescent X-ray emission from SAX 
J1808. 4-3658 was discovered with BeppoSAX on 
1999 March 17-19 (Stella et al. 2000). The ob- 
served MEGS count rate was very low, ^3 x 10^'^ 
cts s~^, which may convert to a luminosity of 
~ 0.7 - 1.5 X 10^2 ergs s"! (0.5-10 kcV) for an 
assumed distance of 2.5 kpc. This is in the range 
of quiescent luminosities in other NS SXTs. 

In this letter, we report results of the 1999 
September ASGA observation of SAX J1808.4- 
3658 made in the quiescent state. 

2. Observation and Results 

SAX J1808.4-3658 was observed with ASGA 
from 1999 September 17, 18:10 (UT) through 
September 20, 4:40, for a net exposure time of 
63 ks. ASGA utilizes grazing-incidence, thin- 
foil mirrors with two kinds of focal plane detec- 
tors: the Gas Imaging Scintillation-proportional 
counter (GIS) and the Solid-state Imaging Spec- 
trometer (Tanaka, Inoue, & Holt 1994, SIS). The 
GIS has a wider field of view, higher time resolu- 
tion, but poorer spatial resolution and low-energy 
efficiency (Ohashi et al. 1996). On the other hand, 
the SIS has a higher spatial resolution, higher en- 
ergy resolution and a better low-energy efficiency. 

For our quiescence observation of SAX J1808.4- 
3658, we had set the GIS in the highest time res- 
olution mode to search for the 401 Hz pulsation 
at the sacrifice of the rise time information and a 
part of the spectral and the spatial information. 



The time resolution was 61 fis in high telemetry 
bit rate and 0.5 ms in medium bit rate. The im- 
age resolution was reduced to 64 x 64 pixels and 
the total number of spectral bins was 256. The 
SIS was operated in 1-ccd faint mode appropriate 
for the observations of dim sources. The SIS data 
have a time resolution of 4 s, and cover a field of 
view of 11' X 11' area. 

Both the SIS and the GIS data were screened 
with the standard data screening criteria. How- 
ever, we could not apply the task 'gisclean', which 
rejects background data using the combination 
of pulse height and rise time information, be- 
cause of the lack of the rise time data. This re- 
sulted in relatively high internal background in 
the GIS data. Using the screened data, we cal- 
culated an image for each sensor. We applied a 
temperature-dependent attitude correction to the 
SIS data, which reduced the systematic error of 
the source position to 0.'5 (Gotthelf 1996"^; Ucda 
et al. 1999b). We found that the images of SIS-0 
and SIS-1 were consistent and both showed slight 
excess fiux at the position consistent with that 
of SAX J1808.4-3658. However, the images of 
GIS-2 and GIS-3 look very different. We show 
in Fig. 1 the images calculated from GIS-2, GIS- 
3 and summed SIS data. The images are slightly 
smoothed with a gaussian kernel appropriate for 
each image. As seen in the GIS images, the dif- 
ference of diffuse emission pattern is apparent, es- 
pecially at the position of SAX J1808.4-3658; a 
point-like emission is seen only in GIS-3. Because 
the image difference (~1 x 10^'^ c/s/cm^/keV) is 
comparable to the local fluctuation of the internal 
background of the GIS (Makishima et al. 1996), 
we tried several different screening criteria using 
the GIS monitor data and the satellite's orbit re- 
lated parameters. However, the difference did not 
disappear. Careful inspection of the GIS-3 image 
showed that the point-like emission at the position 
of SAX J1808. 4-3658 actually consisted of a sin- 
gle bright pixel. Smoothing by a gaussian kernel 
mimicked a point-like source. Because the raw im- 
age does not contain such a bright pixel, the coor- 
dinate transformation to the sky image, which in- 
cludes a position randomization process, may have 
produced such a bright pixel. We confirmed that 
the bright pixel is not produced if we use a differ- 
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ence sequence of random number. Based on these 
analysis, we consider that the coarse image bin- 
ning (1' X 1' bin) and high internal background 
make the GIS data inadequate to search for the 
faint point source. A search for the 401 Hz coher- 
ent pulsations demonstrated also that the statis- 
tics of the GIS data were not enough to detect the 
pulsations, even if the relative pulse modulation 
equals 100 %. 

Because the GIS data proved to be inadequate, 
we used the SIS data to evaluate the X-ray flux 
from SAX J1808.4-3658. We calculated a radial 
profile of the source from the summed SIS image 
(Fig. 2). The center of the profile was taken at 
the local peak of the (smoothed) image, which 
coincided with the position of SAX J1808. 4-3658 
within the accuracy of position determination (0.'5 
in radius). The radial profile was fitted with 
the model radial profile of a point source plus 
a constant background to determine the source 
flux. Precisely speaking, the model radial pro- 
flle slightly depends on the X-ray energy, although 
the dependence is completely negligible compared 
to the statistical errors of the present data. The 
source flux was estimated to be (1.1 ± 0.4) x 10^"^ 
cts s""'^ (average of SIS-0 and SIS-1) for the energy 
range of 0.5-5 keV. Poor statistics of the data did 
not allow us to extract the spectral information. If 
we assume a power law spectrum (photon index of 
2) with an absorption column of 1.3 x 10^^ cm^^, 
which corresponds to the neutral hydrogen column 
of the source direction (Dickey & Lockman 1990), 
then the X-ray flux would be (3 ± 1) x 10"^^ ergs 
s^^ cm^^ in 0.5-5 kcV. This flux converts to a 
luminosity of 1.5 — 3 x 10^^ ergs s^^ for the source 
distance of 2.5 kpc. 

3. Discussion 

We observed the accretion-driven millisecond 
X-ray pulsar SAX J1808.4-3658 with ASCA in 
quiescence. We detected a faint X-ray source at 
the sky position consistent with that of the pulsar. 
However, the source was so dim that we could not 
study its spectral or timing properties. We first 
elaborate on the identification of the source, and 
then on the nature of the source. 

Because the surface density of X-ray sources as 
dim as SAX J1808. 4-3658 is relatively high, an 
X-ray source may happen to coincide with SAX 



J1808. 4-3658 by chance. The chance probability 
may be estimated using the log N— log S relation of 
the X-ray sources. ASCA survey data were used 
to investigate the log N — log S relation for both 
the cosmic X-ray background (Ueda et al. 1999a) 
and the galactic ridge emission (Sugizaki 2000). 
According to those results, extra-galactic sources 
are dominated at a flux level of 3 x lO"^** erg 
s~^ cm~^. The density of sources brighter than 
3 X 10^^"' erg s~^ cm^^ is about 10^ deg^^. Be- 
cause the ASCA position determination error of 
SAX J1808.4-3658 is about 0.5 arcmin, the prob- 
ability to find a source of 3 x 10~^^ erg s~^ cm~^ 
at the position consistent with SAX J1808. 4-3658 
by chance is estimated to be ~2 %. Although this 
probability is not very low, the source we detected 
is most probably SAX J1808.4-3658. If it were 
a chance coincidence of an unrelated source, the 
luminosity we calculated for SAX J1808.4-3658 
should be regarded as an upper limit. 

The quiescent luminosity of SAX J1808.4-3658 
we obtained may be signiflcantly lower than that 
obtained by BeppoSAX (SteUa et al. 2000). In 
general, it is necessary to assume an energy spec- 
trum to convert the flux between different in- 
struments. We assume that the spectral shape 
(eg. photon index, absorption column) of SAX 
J1808. 4-3658 did not change between the ASCA 
and BeppoSAX observations, and only the nor- 
malization did. We have tried two cases for the 
energy spectrum, a power law (P = 2.0) and a 
blackbody (fcT = 0.3 keV), which may be appro- 
priate for NS SXTs in quiescence. Because the 
column density of SAX J1808.4-3658 is not wcU 
constrained, we have tried a range of the column 
density, from 1.3 x 10^^ cm~^ through 6 x 10^^ 
cm^^. Under these assumptions, the MEGS count 
rate (3.1 x lO'^ cts s'^ in 1.3-10 keV) is con- 
verted to the SIS count rate using the utility 
PIMMS, which is developed and maintained by 
HE AS ARC*. Corresponding SIS count rates are 
found to be 0.7 - 1.1 x 10"^ cts s'^ (0.5-5 keV) 
for a power law and 1.9 — 2.9 x 10^^ cts s^^ (0.5-5 
keV) for a blackbody, respectively From this esti- 
mation, we conclude that the observed SIS count 
rate (1.1 ± 0.4 x 10^'^ cts s^^) is at least a factor 
of 4 smaller than the MEGS count rate, even if 
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the statistical errors are considered. This means 
that the X-ray emission from SAX J1808.4-3658 
is time variable. 

The radio pulsar shock emission, which is pref- 
ered by Stella et al. (2000), may explain the time 
variation between the BeppoSAX and the ASCA 
observations. X-ray emission from the shocked 
pulsar wind has been observed from the PSR 
B1259~63/SS 2883 system, which consist of a 48 
ms radio pulsar and a Be star (Hirayama et al. 
1999). ASCA observations of the PSR B1259- 
63/SS 2883 system were made near pcriastron and 
apastron of the highly elliptical orbit (e = 0.86) of 
the pulsar. X-ray emission was detected at both 
periastron and apastron, but the X-ray intensity 
was different by an order of magnitude. Because 
the intensity change is well correlated to the bi- 
nary phase, difference of the stellar wind density 
may cause the intensity change. If the same or 
a similar mechanism works for the X-ray emis- 
sion from SAX J1808. 4-3658, the flux decrease be- 
tween BeppoSAX and ASCA observations may be 
understood as the density decrease of the ambient 
matter around SAX J1808. 4-3658. However, it is 
not clear what would cause the density decrease 
of the ambient matter. Another emission mech- 
anism, neutron star cooling, is also discussed by 
Stella et al. (2000). Although thermal emission 
from the cooling neutron star cannot explain the 
time variation, the mechanism may still be valid 
if the X-ray emission observed by BeppoSAX was 
dominated by another emission mechanism, such 
as residual accretion or radio pulsar shock emis- 
sion. 

Although the BeppoSAX observation favors the 
radio pulsar shock emission for the quiescence 
emission from SAX J1808.4-3658 (Stella et al. 
2000), another possibility can be considered for 
the X-ray emission detected during the ASCA ob- 
servation. ASCA may have observed the X-ray 
emission directly from the radio pulsar. ROSAT 
and ASCA observations of rotation-powered pul- 
sars show that there is a very good correlation 
between the X-ray luminosity and the spin-down 
energy loss rate: Lx = 10~^-Espin (Becker & 
Triimper 1997). Here, the X-ray luminosity, Lx, is 
defined in 0.1-2.4 keV. The observed flux may be 
converted to the unabsorbed luminosity (0.1-2.4 
keV) as 3.5 x lO^^d^ 5 ergs s^^, where ^2.5 is the 
source distance normalized by 2.5 kpc (a power 



law with r = 2 and A^h = 1-3 x 10^^ cm^^ is as- 
sumed). By substituting Lx = 3.5 x lO^^dl 5: we 
obtain the spin-down energy loss rate as -Egpin — 
3.5 X 10^'*d2.5 ^""g^ ^~^- This leads to a surface 
magnetic field of 3.8 x 10*^2.5 G. Although the 
distance to SAX J1808. 4-3658 is not known accu- 
rately, this estimation of the surface magnetic field 
is consistent with the limits obtained from the out- 
burst observation (Wijnands & van der Klis 1998; 
Psaltis & Chakrabarty 1999). This is also consis- 
tent with the surface magnetic field of millisecond 
radio pulsars. Therefore, we consider it possible 
that SAX J1808. 4-3658 might have been a radio 
pulsar during the ASCA observation. 

The very low X-ray luminosity of SAX J1808.4- 
3658 measured by ASCA may not fit in the sce- 
nario that NS SXTs generally have a bottom lumi- 
nosity around 10'^^ — 10^"^ ergs s~^ in quiescence, 
while BH SXTs do not (Asai et al. 1998). We 
argue that the scenario may still be valid, but 
some unique feature of SAX J1808. 4-3658 makes 
it an exception. Because the propeller effect might 
play an important role in NS SXTs in quiescence 
(Menou et al. 1999), one may suspect that the spin 
period or the magnetic field strength of the neu- 
tron star in SAX J1808. 4-3658 have very different 
values than those of other neutron stars in low- 
mass X-ray binaries (LMXBs). However, this is 
not the case. Spin frequencies of the neutron star 
may be estimated from the observations of nearly 
coherent oscillations during type I X-ray bursts 
as observed in some neutron star X-ray binaries 
or from the kHz quasi-periodic oscillations, and 
are mostly found in the range 200-600 Hz (van 
der Klis 2000). According to the recycle scenario 
(Bhattacharaya & van den Heuvel 1991), the sur- 
face magnetic field of the neutron stars in LMXBs 
may be comparable to those of the millisecond ra- 
dio pulsars. Observations of the spectral transi- 
tion of non-pulsating NS SXTs actually gave mag- 
netic field strengths comparable to that of the mil- 
lisecond radio pulsars (Zhang, Yu, & Zhang 1998). 
Therefore, the spin frequency and the magnetic 
field strength of SAX J1808.4-3658 arc not excep- 
tional. It is noteworthy that the outburst luminos- 
ity of SAX J1808.4-3658 was only - lO^^ erg s^^. 
If the quiescent and long-term average outburst 
luminosities are related, as is the case if the qui- 
escent emission is dominated by thermal emission 
from the neutron star surface (Brown, Bildsten, & 



Rutledge 1998; Rutledge et al. 2000), the low qui- 
escent luminosity may result from the low outburst 
luminosity. However, it is not clear why the out- 
burst luminosity of SAX J1808.4-3658 is low (but 
see King 2000). Although we do not know what 
makes SAX J1808. 4-3658 unique, the uniqueness 
may be responsible for both the detectable spin 
modulation during the outburst and the very low 
luminosities in quiescence. 
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Fig. 1. — X-ray images calculated from the ASCA data; (a) GIS-2, (b) GIS-3, and (c) sum of SIS-0 and 
SIS-1. Expected position of SAX J1808. 4-3658 is indicated by the cross with a circle. GIS images are 
smoothed with a gaussian kernel of cr = 1'. Contours are plotted in linear scale from 13.7 cts arcmin"^ with 
an increment of 1.37 cts arcmin^^. SIS image (sum of SIS-0 and SIS-1) is smoothed with a gaussian kernel 
(cr ~ 0.42') after rebinning by a factor of 8. Contours are plotted in linear scale from 63.3 cts arcmin^^ with 
an increment of 4.87 cts arcmin"^. Note that the field of view is different between GIS and SIS. 
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Fig. 2. — Radial profile of the X-ray image calculated from the sum of the SIS-0 and SIS-1 data. The local 
peak of the image was taken as the center of the radial profile. Data points are plotted as crosses and the 
best-fitting model profile for a point source (plus a constant background) as a sold line. No background is 
subtracted from the data. 



